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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  NOTE  NO.  1854 

APPROXIMATE  CORRECTIONS  FOR  THE  EFFECTS  OF  COMPRESSIBILITY 
ON  THE  SUBSONIC  STABILITY  DERIVATIVES  OF  SWEPT  WINGS 
By  Levis  E.  Fisher 

SUMMARY 


A  method  which  has  previously  heen  presented  for  the  estimation  of 
the  low- speed  stability  derivatives  of  swept  wings  is  modified  by  means 
of  the  Prandtl-Glauert  rule  to  yield  approximate  corrections  for  the 
first-order  three-dimensional  effects  of  compressible  flow  in  the  sub¬ 
sonic  region.  Corrections  are  presented  in  chart  form  for  several  aspect 
ratios  and  angles  of  sweep  and  may  be  computed  rather  simply  for  other 
aspect  ratios  and  angles  of  sweep. 

Lift-curve  slopes  and  values  of  damping  in  roll  are  compared  with 
those  obtained  by  use  of  a  generally  accepted  method  and  are  discussed  in 
relation  to  unpublished  data. 


INTRODUCTION 


The  three-dimensional  effects  of  compressibility  on  the  lift-curve 
slopes  of  unswept  wings  are  treated  in  references  1  and  2,  and  the  treat¬ 
ment  is  extended  to  the  slopes  of  swept  wings  in  references  3  and  4.  The 
only  other  stability  derivative  that  appears  to  have  been  considered  is 
the  damping  in  roll,  for  which  compressible-flow  values  can  be  estimated 
from  reference  5  • 

At  present  no  rigorous  theories  exist  that  seem  suitable  for  the 
evaluation  of  all  the  stability  derivatives  of  swept  wings  even  in  incom¬ 
pressible  flow.  The  effects  of  compressibility  on  these  derivatives  may 
not  be  treated  in  an  exact  manner  for  some  time.  For  many  applications, 
however,  such  as  in  dynamic  stability  calculations,  an  approximation  of 
the  effects  of  compressibility  is  adequate. 

A  complete  analysis,  based  on  strip  theory  modified  by  simple  means 
to  account  for  the  primary  effects  of  aerodynamic  induction,  has  been 
presented  for  the  contribution  of  sweep  to  the  stability  derivatives  of 
wings  for  zero  Mach  number.  (See  reference  6.) 
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The  method  forms  a  convenient  basis  for  making  estimates,  in  a 
rational  manner,  of  compressibility  effects  on  the  stability  derivatives. 

This  estimate  is  accomplished,  in  the  present  paper,  by  increasing  the 
lift-curve  slope  of  each  wing  section  by  an  application  of  the  Prandtl- 
Glauert  rule  and  by  using  the  equations  of  reference  6  to  derive  correc¬ 
tions  for  the  first-order  three-dimensional  effects  of  compressibility 
in  the  subsonic  Mach  number  region.  The  corrections  are  applied  in  the 
form  of  ratios  of  the  compressible-flow  equations  to  the  incompressible- 
flow  equations;  thus,  the  effects  of  errors  in  the  theory  are  minimized. 

The  same  procedure  is  applied  to  all  the  steady- state  derivatives 
with  the  exception  of  the  yawing  moment  due  to  yawing  and  the  lateral 
force  due  to  yawing.  Since  the  contribution  of  a  wing  alone  to  either 
of  these  derivatives  is  small  compared  to  that  of  the  vertical  tail,  Mach 
number  effects  on  these  derivatives  would  appear  to  be  insignificant. 

Consideration  of  the  lift-curve  slope  and  the  damping  in  roll  is 
included  herein  for  completeness,  although  more  rigorous  procedures  for 
treating  these  derivatives  are  discussed  in  references  1  to  5 • 


SYMBOLS 


The  symbols  used  in  the  analysis  and  in  the  presentation  of  results 
are  defined  herein.  All  spans  and  chords  are  measured  perpendicular  and 
parallel,  respectively,  to  the  plane  of  symmetry. 


lift  coefficient  /  LlftN 

|pv2s 


rolling-moment  coefficient  [Bolling  moment 


|pV2Sh 


m 


pitching-moment  coefficient  /^Itching  moment 


fpV2Sc 


n 


yawing-moment  coefficient  /Yawing  moment 


pV2Sh 


Cy 


lateral-force  coefficient  /Lateral  force 
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A 

A 

M 

A 

a 


section  primary-force  coefficient 
( see  reference  6) 

mass  density,  slugs  per  cubic  foot 
wing  area,  square  feet 
free- stream  velocity,  feet  per  second 
wing  span,  feet 
wing  chord,  feet 
wing  mean  chord,  feet  (S/h) 

longitudinal  distance  rearward  from  airplane  center  of  gravity 
to  wing  aerodynamic  center,  feet 

spanwise  distance  from  plane  of  symmetry  to  any  station  on 
wing  quarter- chord  line,  feet 

spanwise  distance  from  plane  of  symmetry  to  effective  lateral 
center-of -pressure  location  of  resulting  load  causing 
rolling  moment,  feet 

aspect  ratio  (b2/s) 

taper  ratio  /Tip  chord  \ 

\R00t  chord/ 

Mach  number  f _ 

\Speed 

angle  of  sweep  measured  at  quarter-chord  line,  degrees 
angle  of  attack,  radians 
angle  of  sideslip,  radians 

local  angle  of  sideslip;  angle  between  plane  of  symmetry  and 
local  air- stream  direction  at  quarter-chord  point  of  any 
section,  radians 


1 _ ) 

of  sound/ 


Section  primary  force' 
|pV2c 
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section  lift-curve  slope  for  section  normal  to  quarter- chord 
line  when  placed  in  direction  of  free  stream 

rolling  angular  velocity,  radians  per  second 

pitching  angular  velocity,  radians  per  second 

yawing  angular  velocity,  radians  per  second 


£b 

27 


nondimensional  rolling- velocity  parameter 


22 

27 

rb 

27 


nondimensional  pitching-velocity  parameter 
nondimensional  yawing-velocity  parameter 


Subscripts : 


L 

E 

M 

A=0° 


a,p,p,r,q 


left  wing  panel 
right  wing  panel 
value  at  Mach  number  M 
restricted  to  zero  sweep 

denotes  derivative  of  CL,  Cm,  Cj,  Cn,  or  Cy  with  respect 
t0  3’  W  27;  f0r  e3Cample^  %  B 


ANALYSIS 


In  reference  1,  three  variations  of  the  application  of  the  theory  of 
small  perturbations  to  determine  the  effects  of  compressibility  are  sum¬ 
marized.  Two  of  these  methods  consist  of  changing  the  dimensions  of  the 
wing  or  the  angle  of  attach  or  both  through  the  Prandtl-Glauert  trans¬ 
formation  and  then  obtaining  the  aerodynamic  characteristics  of  the 
resulting  wing  from  incompressible-flow  theory .  These  concepts  are  dis¬ 
advantageous  in  that  the  plan  form,  of  the  wing  in  sideslip  or  yawing  iB 
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distorted  unsymmetrically •  The  third  method,  which  consists  of  increasing 
the  two-dimensional  lift-curve  slope  of  a  wing  for  which  the  dimensions 
were  held  constant,  is  the  m®st  convenient  to  use  and  is  the  procedure 
adopted  herein . 

Eeference  1  shows  that  the  increase  in  section  lift-curve  slope  is 
the  only  first-order  three-dimensional  effect  of  compressibility  on  a  thin 
wing.  For  an  unswept  and  untapered  wing  in  a  free- stream  flow  of  Mach 
number  M,  the  section  lift-curve  slope  aQ  is  shown  to  he  increased  by 
the  Prandtl-Glauert  factor  l/B0  where  B0  =  n/T  -  .  If  the  wing  is 

swept  or  oblique  to  the  flow,  however,  only  the  component  of  the  Mach 
number  in  the  direction  normal  to  the  quarter- chord  line  is  effective  in 
increasing  the  section  lift-curve  slope  of  the  wing.  (See  reference  7.) 

Thus  the  lift  of  a  thin  two-dimensional  wing  which  is* oblique  to  the  flow 

by  an  angle  A  is  increased  by  the  factor  l/B  where  B  =  \/l  -  M2cos2A  • 

A  somewhat  greater  rate  of  increase  of  wing  lift  with  Mach  number  would  be 
expected  for  wings  of  finite  thickness  (reference  8),  but  it  is  shown  in 
reference  k  that  this  effect  is  small  for  usual  aspect  ratios  and  sweep 
angles.  For  the  present  analysis,  therefore,  the  effects  of  thickness 
will  be  neglected. 

The  stability  derivatives  under  consideration  in  this  analysis  depend 
substantially  upon  the  orientation  and  the  rate  of  change  of  lift  and 
induced-drag  (which  depends  on  lift)  forces  due  to  a  change  in  attitude  of 
the  wing.  When  these  lift  forces  are  increased  in  accordance  with  the 
Prandtl-Glauert  rule,  the  stability  derivatives  that  depend  so  directly 
upon  the  lift  forces  may  be  expected  to  include  the  first-order  effects  of 
compressibility.  The  equations  of  reference  6  consider  the  two-dimensional 
effects  of  sweep  corrected  for  finite  aspect  ratios.  If  the  two-dimensional 
effects  are  further  corrected  for  compressibility,  the  resulting  equations 
should  also  be  valid  for  three-dimensional  flow. 

In  reference  6,  strip-theory  equations  were  first  derived  for  each  of 
the  derivatives  with  approximate  consideration  given  to  aerodynamic  induc¬ 
tion.  Since  the  strip  theory  would  be  quantitatively  in  error,  these 
equations  were  used  only  to  obtain  corrections  to  unswept-wing  derivatives, 
which  could  be  calculated  by  more  exact  theories . 

In  the  derivations  presented  in  reference  6,  the  section  lift-curve 
slope  was  left  arbitrary  only  when  it  had  a  primary  effect  on  the  deriva¬ 
tive;  that  is,  when  the  derivative  was  proportional  to  aQ .  In  other 

instances  aQ  was  assumed  equal  to  2tt.  For  the  determination  of  the 
compressibility-correction  equations  in  the  present  paper,  the  strip- 
theory  equations  of  reference  6  are  restated  in  more  basic  form  by 
retaining  a0  in  all  instances.  The  compressible-flow  equation  is  estab¬ 
lished  by  substituting  a0  in  the  incompressible-flow  development  with  a0/B 
and  then  by  replacing  aQ  with  2m.  In  the  case  of  sideslipping  or  yawing, 

the  problem  seems  a  little  more  complex  unless  it  is  understood  that  the  local 
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conditions  at  each  wing  section  are  considered.  The  Mach  number  correction 
is  then  the  ratio  of  the  compressible -flow  equation  to  its  incompressible - 
flow  counterpart.  This  ratio,  as  it  variea  with  Mach  number,  may  readily 
be  applied  as  a  correction  to  incompressible -flow  values  of  a  derivative 
which  may  be  determined  from  reference  6  or  from  experimental  data.  The 
derivations  for  the  Mach  number  corrections  for  all  stability  derivatives 
considered  are  given  in  the  appendix.  Compressibility  effects  on  the 
derivatives  and  Cy^,  are  not  considered  in  this  paper,  although  a 

logical  method  for  the  estimation  of  the  correction  to  Cn  is  suggested 

in  the  appendix. 

The  present  method  is  subject  to  the  limitations  of  the  linear  pertur¬ 
bation  theory  that  assumes  only  small  departures  of  the  fluid  velocity 
from  the  free-stream  velocity.  The  strip  theory  of  reference  6  is  limited 
for  most  accurate  results  to  high  aspect  ratios  and  taper  ratios  close  to 
unity.  Errors  inherent  in  the  equations,  ^such  as  those  magnified  by  low 
aspect  ratios,  however,  tend  to  nullify  each  other  when  the  correction  is 
applied  in  the  form  of  a  ratio.  For  the  same  reason  the  equations  are 
considered  to  be  applicable  to  wings  that  are  tapered  at  least  moderately . 


RESULTS  AND  DISCUSSION 


The  correction  factors  derived  for  some  of  the  derivatives  depend  on 
the  static  margin  x/c-  In  order  to  provide  an  indication  of  the  impor¬ 
tance  of  the  static  margin,  the  factors  dependent  upon  it  were  calculated 
for  static  margins  of  0  and  0.2  for  wings  having  sweep  angles  of  30° 
and  50°  and  an  aspect  ratio  of  4 .  The  results  shown  in  figure  1  indicate 
that  the  static  margin  is  of  little  importance;  consequently,  the  margin 
was  assumed  to  be  zero  for  those  equations  in  which  it  appeared. 

The  corrections  which  may  be  multiplied  by  the  incompressible-flow 
values  of  the  stability  derivatives  under  consideration  in  order  to  arrive 
at  approximate  values  for  the  derivatives  at  Mach  numbers  between  zero  and 
unity  are  presented  in  chart  form  (figs.  2  to  9) •  The  charts  present  cor¬ 
rections  for  four  aspect  ratios  (2,  3>  4,  and  6)  and  five  angles  of  sweep 
(0°,  30°,  40°,  50°,  and  60°)  .  Compressibility  corrections  for  aspect 
ratios  and  sweep  angles  other  than  for  those  presented  may  be  obtained 
either  by  interpolation  of  the  given  curves  or  by  calculation  from  the 
given  equations,  (See  table  I .)  These  corrections  are  applicable  to  either 
sweptback  or  sweptforward  wings  within  the  limitations  of  the  basic  theory . 

The  calculated  effects  of  compressibility  on  the  stability  derivatives 
of  representative  wings  having  an  aspect  ratio  of  4  are  presented  in 
figure  10,  which  shows  that  an  Increase  in  sweep  may  be  expected  to 
decrease  the  magnitude  of  the  increments  due  to  compressibility.  Values  of 
the  derivatives  for  zero  Mach  number  were  obtained  from  reference  6  and  any 
effects  of  static  margin  were  neglected. 
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A  comparison  of  the  present  method  with  the  generally  accepted  method 
of  equivalent  wings  for  determining  approximate  Mach  number  effects 
on  Or  and  C7  is  given  in  figures  11  to  14.  The  use  of  the  equi  valent - 

6p 

wing  method  in  conjunction  with  calculations  made  hy  the  Weissinger  method  is 
suggested  in  references  3  and  5  for  CL  and  ,  respectively.  A  pre¬ 
liminary  comparison  of  these  methods  with  some  unpublished  experimental 
data  indicates  that  both  methods  underestimate  the  effects  of  compressi¬ 
bility  on  Cq^,  particularly  for  the  smaller  sweep  angles  and  higher  Mach 

numbers.  This  underestimation  is  probably  due  in  part  to  the  neglect  of 
the  finite -thickness  effect  of  the  experimental  wings.  The  present  method 
indicates  greater  variations  with  Mach  number  than  the  equivalent-wing 
method  (see  figs.  11  to  Ik)  and,  therefore,  is  expected  to  give  better 
agreement  with  experimental  results . 

The  present  method  was  derived  by  application  of  a  form  of  strip 
theory  to  wings  having  a  taper  ratio  of  1.0.  An  indication  of  the  reli¬ 
ability  of  the  method  for  moderately  tapered  wings  is  provided  for  the 
derivatives  and  Cj  in  figures  11  and  13,  respectively.  These 

figures  show  that  for  an  aspect  ratio  of  4  the  differences  between  the  pre¬ 
sent  method  and  the  method  of  equivalent  wings  are  about  the  same  for  taper 
ratios  of  either  1.0  or  0.5*  The  present  method,  therefore,  might  be 
expected  to  apply  with  reasonable  accuracy  to  moderately  tapered  wings . 
Comparisons  of  the  present  method  with  the  equivalent-wing  method  for 
untapered  wings  having  aspect  ratios  of  2  and  6  show  that,  over  this  range 
of  aspect  ratios,  the  present  method  yields  consistently  larger  values  of 

the  corrections  for  either  Ct  or  C7  .  (See  figs.  12  and  14.)  Both 

La  "p 

methods  are  in  general  agreement  in  that,  for  a  given  aspect  ratio,  smaller 
effects  of  compressibility  are  indicated  as  the  sweep  angle  is  increased. 


CONCLUDING-  REMARKS 


An  adaptation  of  the  Prandtl-Glauert  rule  is  used  to  modify  existing 
theory  for  the  first-order  effects  of  compressibility  on  the  subsonic  sta¬ 
bility  derivatives  of  swept  wings.  Lift-curve  slopes  and  values  of  damping 
in  roll  are  compared  with  those  obtained  by  use  of  a  generally  accepted 
method  and  are  discussed  in  relation  to  unpublished  data.  Because  of  the 
lack  of  experimental  data,  the  reliability  of  the  present  method  for  the 
calculation  of  the  effect  of  high  subsonic  Mach  numbers  on  all  the  sta¬ 
bility  derivatives  is  not  known ^  however,  the  method  is  expected  at  least 
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to  provide  a  reliable  indication  of  the  order  of  magnitude  of  the  Mach 
number  effect. 

In  general,  the  effects  of  compressibility  appear  to  become  smaller 
as  the  angle  of  sweep  is  increased . 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 

Langley  Air  Force  Base,  Va.,  January  20,  19^9 
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APPENDIX 

DERIVATION  OF  COMPRESSIBILITY-CORRECTION  EQUATIONS 
Lift-Curve  Slope 


The  strip-theory  relationship  for  fC T  'N  in  its  "basic  form,  as 

„  V  a/M=0  ’ 


given  "by  reference  6,  is 


aQ  cos  A 
aQ  cos  A 


The  application  of  the  Prandtl-Glauert  rule  gives  the  result 


(Clo)m  1 


-  cos  A 

ao  cos  A 
B  7rA 


If  the  thin -  airfoil  -  the  ory  value  of  27T  is  now  assigned  to  a  ,  the 
compressibility-correction  equation  becomes 


(Clo)m=i 


A  +  2  cos  A 
AB  +  2  cos  A 


At  zero  sweep  angle,  equation  (l)  reduces  to  the  equation  derived  in 
reference  2  by  increasing  the  lift  coefficient  hy  i/b2  and  decreasing 
the  span  and  angle  of  attack  "by  B  in  the  lifting-line-theory  equation 
for  CT  . 

Equation  (l)  indicates  an  increase  in  the  compressibility  factor 
from  unity  to  l/B  as  the  aspect  ratio  increases  from  zero  to  infinity. 

Foiling  Flight 

Foiling  moment  due  to  rolling.-  A  procedure  similar  to  that  employed 
for  CL^  can  "be  used  to  develop  an  expression  for  the  effects  of  compres¬ 
sibility  on  the  damping  in  roll  based  on  the  following  equation  from 

reference  6  in  basic  form: 


(C?p)m=( 


Aa0  cos  A 
^  2a0  cos  a 
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The  application  of  the  Prandtl-Glauert  rule  gives 

A  —  cos  A 

fn  \  —  _  -l-  _ 

^  zpjM  ”  Z  a  +  £0  2  cos  A 
B  T 

If  aG  =  27T,  then  the  compressibility-correction  equation  is 


V^p/M  _  A  +  b  cos  A  (2) 

(%)M=0  **  +  4  008  A 

Yawing  moment  due  to  rolling.-  The  incompressible-flow  relationship 
for  C^,  based  on  reference  6,  is 

„  •  .  .  tL  .  aD  cos  A \fi  tan  A  +  tan2A\ 

.  CL  A  *  HA  *  — 5r^ A*  -A"  —1S7 

V  npjM=0  -  6  2a.  00s  A 

A  +  — 2 - 

7T 


When  modified  for  compressibility,  the  equation  becomes 


(%)m  =  "  T 


c(.  _l.  a°  008  a\(*  tan  A  +  tan-2A ) 

V  ~B~aF~A°~A~  12  / 


2  cos  A  O 


By  letting  aQ  =  2rr,  the  compressibility-correction  equation 


becomes 


— ) 

cLM=0 


A  +  b  cos  A 
AB  +  ^  cos  A 


AB  +  6[AB  +  cos  A 


tan  A 


tan2/y\ 
12  / 


A  +  6\A  +  cos  aJ(^5 


'  x  tan  A  .  tan: 


Fieure  1  shows  that  the  static  margin  x/c  has  a  negligible  effect 
due  to  rolling  becomes 
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AB  +  |UB  +  cos  Ajtan2A 
A  +  4  cos  A  _ j _ {_ _ 

AB  +  4  cos  A  1  /  \  o 

A  +  —  fA  +  cos  AJtanM 


(3) 


Lateral  force  due  to  rolling The  equation  for  Cy  ,  as  developed 

IP 

for  the  incompressible -flow  analysis  of  reference  6,  is 


(%)m-o 


aD  cos  A 
A  +  - 

=  Ct  tan  A - — - 

A  +  2ao  coa  A 
7 r 


The  value  for  the  derivative  at  Mach  number  M  may  be  written 


A  + 


a 


Cp  tan  A 


_o  cos  A 

B  27 r 


A  + 


'O  2  cos  A 


B 


7T 


If  a0  =  27 r,  then 


A  +  k  cos  A  AB  +  cos  A 
AB  +  4  cos  A  A  +  cos  A 


Sideslipping  Flight 

Bolling  moment  due  to  sideslip.-  In  the  case  of  a  swept  wing  in 
sideslip,  the  different  effective  sweep  angles  produce  a  difference  in 
effective  Mach  number  between  the  left  and  the  right  wing  panels .  Cal¬ 
culations,  however,  show  this  difference  in  Mach  number  to  be  of  second- 
order  importance .  The  effect  may,  therefore,  be  neglected  and  a  proce¬ 
dure  similar  to  that  employed  heretofore  may  again  be  used. 


The  equation  for 
reference  6,  is 


as  developed  for  incompressible  flow  in 


12 
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( 


CL/M=0 


A  + 


a0  cos  A 
7 r 


JJB  ]  _ _ 7 T  tan  A 

CL/a=0°  */2  A  +  gao  COB  A  2 


(5) 


TT 


The  procedure  employed  In  the  present  paper  is  actually  applicable 

only  to  the  increment  in  that  is  due  to  sweep.  Although  the  other 

increment  cannot  he  handled  by  the  present  simplified  proced^e  itprob- 
ably  is  affected  hy  Mach  number  in  much  the  same  manner  as  the  incremen 
due  to  sweep,  since  both  increments  result  from  load 

same  general  form.  The  compressibility  factor  determined  f or the  incre 
ment  due  to  sweep,  therefore,  is  assumed  to  apply  to  the  total  value  of 
the  derivative  .  The  additional  assumption  is  made  that  compressibility 
causes  no  material  change  in  the  lateral  center  of  pressure . 


If 


(%)  iS  neglected  and  the  Prandtl-Glauert  rule  is  applied 

V  cL/a=0° 


to  the  increment  due  to  sweep,  the  equation  becomes 


,  v  A  +  — 

/C7  \  -  B 

("c^L  - "  $2  775 


ao  cos  A 


7 r 


®o  2  cos 
B 


_  tan  A 
A  2 


TT 


If  the  thin-airfoil-theory  value  of  2tt  is  substituted  for  aQ,  the 
compressibility  factor  becomes 


A  +  4  cos  A  AB  +  2  cos  A 
AB  +  4  cos  A  A  +  2  cos  A 


(6) 


Equation  (6)  should  now  be  employed  in  the  following  manner  to  cor 
rect  (Cl$M=0  as  calculated  from  equation  (5); 

~  co£  A 


A  + 


ClA=0°  b^2 


7 r 


tan  A 


A  + 


2a0  cos  A  2 


7 T 


A  +  4  cos  A  AB  +  2  cos  A 


AB 


+  k  cos  A  A  +  2  cos  A 
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where  ( _ §.]  and.  .7  may  "be  determined,  from  reference  6 

V  Cl/A=0°  W2 


Yawing  moment  due  to  sideslip.-  The  equation  for  as 

developed  in  reference  6,  is 


Pl2/m=o 


cl2/a=o° 


tan  a  /x  ao  sln  A  +  ao  cos  A 

/  2a0  cosA\\a  TrA  2ir 

7TA(A  +  — 2 - 


A  _  ArTT 
2  4a0  cos  A 


If  — £  is  neglected  and  the  Prandtl-Glauert  rule  is  applied 

\  L /A=0° 

to  the  increment  due  to  sweep,  the  equation  "becomes 


'L  /  M 


_ tan  A 

TrA  (a  +  f2  2  cos  A 
V  B  7T  , 


i  sin  A  +  fo  cos  A A 
c  ttA  B  2tt  2 


cos  A 


By  letting  aQ  =  2ir,  the  compressihility  factor  is 


,CL2/M 


cl2/m=o 


A  +  4  cos  A  c  A 
AB  +  4  cos  A 


ci  sin  A  ,  „  AB 

0=-  -  +  cos  A — 


cy.  sin  A  .  „„„  .  A 
%  “A-  +  008  A  -  - 


A2B2 
8  cos  A 


A 

cos  A 


If  x/c  is  restricted  to  zero,  the  equation  "becomes 0. 


v 
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CL7M 


CL7M=0 


A  +  4-  cos  A  aS2  +  4AB  cos  A  -  8  cob2A 

AB  +  4  cos  A  .p.i.  o  ?a 

Ad  +  4A  cos  a  —  8  coswV 


Equation  (8)  may  now  be  used  with  respect  to 


_ ]  in  the 


cl2/m=o 


manner  as  equation  (6)  was  used  with  respect  to  ) 

\  CL /  M=0 


The  deriv¬ 


ative 


may  he  determined  from  reference  9 


L  /A-0 


Lateral  force  due  to  sideslip.-  The  basic  equation  for 
developed  in  reference  6,  is 


h)  _  3a0  tan  A  sin  a 

CWM=0  ^  /  +  2ap  cos  Av 
\  TT 


The  corresponding  Mach  number  equation  becomes 


tan  A  sin  A 


cl2/m 


TT^AlA  +  —  2_cos_A 


B  7 T 


By  letting  aQ  =  27T,  the  compressibility  factor  is 


cl2/m 


A  +  4  cos  A 
AB  +  4  cos  A 


CL  /M=0 
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V 


If  x/c  is  restricted  to  zero 


A^tan^A  +  2. 
AB  +  6  cos  A  B 

A^tan^A  +  ^ 
A  +6  cos  A 


(11) 


Yawing  Flight 

E oiling  moment  due  to  yawing.-  In  the  determination  of  the  effect  of 
compressibility  on  the  yawing  derivatives  of  a  wing,  another  factor  is 
encountered  which  must  be  considered;  that  is,  the  spanwlse  variation  in 
Mach  number. 

In  reference  6,  Ci  1®  shown  to  be  the  result  of  two  force  coef¬ 
ficients  c1l  and  c^  acting  on  the  left  and  right  wings,  respectively . 

The  following  equation  is  presented: 


w 
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I 


k 


« 


In  equations  (13),  the  expression  In  the  first  parenthesis  is  related  to 
the  wing  lift-curve  slope,  which  depends  only  on  the  free-stream  Mach 
number;  the  expression  in  the  second  parenthesis  is  related  to  the  local 
wing  section  lift  coefficient,  which  varies  along  the  span  and  must  he 
modified  for  the  spanwise  variation  of  Mach  number;  the  expression  in 
the  third  parenthesis  modifies  equations  (13)  for  the  spanwise  variation 
in  velocity  in  yawing  flow. 


By  letting 


B, 


=  'y/l  -  M2cos2A^1  +  ~r 


~  B- 


M^cos^A 

B 


and 


Bp  -  y/l  -  M2cos2A^1  —  2^ 


Sf  B  + 


M^cos^A 

B  ~ 


equations  (13)  -when  modified  for  compressibility  become 


cos  A 


p »  -2  sin  A 


1  + 


o  cos  A 
Bp  7rA 


a, 


1  +  2 


o  cos  A 


Bp  ttA 


and. 


c*i  _  C. 
"  - 


1+^=2  cos  A 

B  TTA 


B 


cos  A 


—  cos  A 

Vr 


ao 

P  1  sin  A 


"R  + .  L 

-1  +  £os  ,.A  ,  „ao  cos  A  I  V1 

}+3R  TTA  1  +  \  ~Wa~I 


where,  from  reference  6 


p’  =~§  [b72 +  (b^-i)6 


an  A 


* 
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+  tan  A 


ac  cos  A 

7rA  5  (  J 

2an  cos  A  b/2  lb/2 


tan  A 


After  substituting  this  value  of  c-j_  —  C]_  from  equation  (lb)  into 

equation  (12),  integrating  across  the  span,  differentiating  the  result 
vith  respect  to  the  yaving- velocity  parameter  and  replacing  aQ 

by  27 r,  the  Mach  number  value  of  the  derivative  of  rolling  moment  due  to 
yaving  becomes 


1±\  -  I  o  A(1  -  B2)  AB  +  2  cos  A  /x  tan  A  tan^A 

ClL  “  6  d  +  B(AB  +  2  cob  A)  AB  +  %  cos  A  l  c  2A  24, 


At  zero  Mach  number  this  equation  reduces  to  the  equation  presented  in 
reference  6: 


^r  |  _  1  A  +  2  cos  A  /x  tan  A  +  tangA 
clA=0  3  A  +  k  cos  A  \c  2A  2k , 
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Yawing  moment  due  to 


and  lateral  force  due  to  yawing The 


contribution  of  a  wing  alone  to  either 


Is  small  compared 


to  the  contribution  of  the  vertical  tailo  For  this  reason,  these  two 
derivatives  were  not  deemed  as  being  of  sufficient  importance  to  warrant 
consideration  similar  to  that  given  the  derivative  The  derivatives 


may  be  corrected  for  compressibility  effects  on  a  similar  basis;  however, 
the  procedure  used  for  these  derivatives  could  not  be  expected  to  be  very 
reliable  because  the  drag  and  the  drag  distribution  are  important  but  can¬ 
not  be  handled  in  a  logical  manner »  Perhaps  the  most  reliable  procedure 
for  the  derivative  Cn  is  the  use  of  the  incompressible- flow  equation 

of  reference  6  with  the  wing  profile-drag  coefficient  appropriate  to  the 
Mach  number  in  question  ■ 
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,ch  number  for  a  representative  wing  with  two  angles  of  sweephack. 
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Corrections  for  the  effect  of  compressibility  on  the  lift— curve  slope^,  equation  (l) 

the  lift  due  to  pitching ,  equation  (10) . 
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igure  3.—  Corrections  for  the  effect  of  compressibility  on  the  rolling  moment  due  to  rolling 

equation  (2),  and  the  lateral  force  due  to  sideslip^  equation  (9). 
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Figure  km—  Corrections  for  the  effect  of  compressibility  on  the  yawing  moment  due  to  rolling 

equation  (3).  =  0. 
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Correct  ions  for  the  effect  of  compressibility  on  the  lateral  force  due  to  rolling, 

equation  (k) . 


(c/pg) 


NACA  TN  No.  1854 


Corrections  for  the  effect  of  compressibility  on  the  rolling  moment  due  to  sideslip, 

equation  (6). 
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Corrections  for  the  effect  of  compressibility  on  the  yawing  moment  due  to  sideslip 


32 


NACA  TN  No.  I85U 


\ 

\\ 

V 

( 

\ 

\ 

\ 

§  §  § 


N 

N  ' 

V 

\\ 

1 

\ 

[\ 

\ 

0 

\ 

' 

u 

^ 


0 

1 

°o 

it 

A 

\ 

\ 

j 

1 

I r 

□ 

□ 

S’ 

1 1 

I  * 

^  'b 

Qj  £ 

<0  ^ 

s;  I 

q:  ^ 


§§  §  3 


=rpc 

\v 

=p! 

-§-P-  -8 

tl 

I  11 1  u  x 

cr\  \o  v°  _ 
«  «  <3  ® 


5 

0 

<vb 

\\ 

I, 

"V 

V 

r 

r 

\ 

1 

b_ 

V 

11 

§  ^  § 


§  8 


'■Q  r~ 


^iS<s 


r— I 

O 

rf 

,3 

O 

-P 

0 

03 

a 

•P 

O 

0 

P 

ch 

ch 

rd 

0 

0 

-P 

U 

cd 

0 

& 

d 

0 

R 

1— l 

PJ 

ed 

0 

■d 

O 

CQ 

0 

d 

1 — 1 

W 

cd 

ri 

J> 

•H 

-P 

> 

cd 

0 

3 

r_j 

V 

0 

I 

i — i 

0 

od 

H 

0 

rQ 

•H 

P 

0 

0 

0 

ch 

0 

P 

CQ 

Qi 

»d 

& 

O 

O 

rR 

O 

"P 

O  u  00 
o 

.  u  © 

I  0  o 
•  Rh  ri 

H  cd  0 
1 — I  -P  P 
0 

0  O  Ch 
h  >  P 

d  -P  !h 


* 


Figure  12.—  Comparison  of  two  methods  for  calculating  Mach  number  effects  on  lift— curve  slopes  for 
two  aspect  ratios,  X  =  1.0.  (incompressible— flow  values  calculated  by  method  of  Weissinger  in 
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ratios.  A  =  4.  (incompressible—: flow  values  calculated  by  method  of  Weissinger  in 
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